Onset of collective and cohesive motion 
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We study the onset of collective motion, with and without cohesion, of groups of noisy self- 
propelled particles interacting locally. We find that this phase transition, in two space dimensions, 
is always discontinuous, including for the minimal model of Vicsek et al. [Phys. Rev. Lett. 75,1226 
(1995)] for which a non-trivial critical point was previously advocated. We also show that cohesion 
is always lost near onset, as a result of the interplay of density, velocity, and shape fluctuations. 
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Collective motion can be observed at almost every scale 
in nature, from the familiar human crowdsQ, bird flocks 
and fish schools 2], to unicellular organisms like amoe- 
bae and bacteria |3(, individual cells and even at mi- 
croscopic level in the dynamics of actin and tubulin fil- 
aments and molecular motors [a, |fj . Whereas biologists 
tend to build detailed representations of a particular case, 
the ubiquity of the phenomenon suggests underlying uni- 
versal features and thus gives weight to the bottom-up 
modeling approach usually favored by physicists Q. 

In this respect, the simple model introduced by Vicsek 
and collaborators stands out because of its minimal 
character and the a priori least-favorable conditions in 
which it is defined. In the Vicsek model (VM), identical 
pointwise particles move at constant velocity and inter- 
act locally by trying to align their direction with that of 
neighbors. Remarkably, even in the presence of noise and 
in the absence of leaders and global forces, orientational 
long-range order arises, i.e. collective motion emerges, if 
the density of particles is high enough or, equivalently, if 
the noise is weak enough. The existence of the ordered 
phase was later "proved" by a renormalisation-group ap- 
proach based on a phenomenological mesoscopic equa- 
tion [9j. More recently, this work was extended to the 
case where the ambient fluid is taken into full account, 
yielding novel mesoscopic equations for suspensions of 
self-propelled particles [lflf. 

The nature of the non-equilibrium phase transition 
to collective motion, however, is not well established. 
Vicsek et al. concluded from numerical simulations in 
two and three dimensions that it is continuous ( "second- 
order") and characterized by a set of critical indices, but 
these results remain somewhat crude ^]], even though 
the undeniably minimal character of the VM makes it a 
good candidate for representing a universality class. 

Moreover, from a modeling point of view, an often de- 
sirable ingredient missing in the VM is cohesion: when 
put together in an infinite space, particles do not stay to- 
gether and fly apart. In other words, no collective motion 
is possible in the zero-density limit of the VM. Recently, 
we have shown how one can ensure cohesion in simple 
models derived from the VM without resorting to leader 



particles or long-range or global forces [l2j |. 

In this Letter, we study the onset of collective motion 
with and without cohesion in this very general setting, 
trying to assess the universality of the results of Vicsek 
et al. In both cases, we find that the onset of collec- 
tive motion in the VM and related models is actually 
discontinuous ("first-order") and that its apparent con- 
tinuous character is due to strong finite-size effects. We 
also show that without cohesion, the transition point is 
nevertheless accompanied by a non-trivial superdiffusive 
behavior of particles which, we argue, could be measured 
experimentally. In the presence of cohesion, our study 
reveals that the onset of collective motion is the theater 
of a complex interplay between density, velocity, sound 
and shape modes, giving rise to fascinating dynamics. 

The original VM is defined as follows: identical point- 
wise particles move synchronously at discrete timesteps 
At — 1 by a fixed distance vq. In two space dimen- 
sions — to which we restrict ourselves in the following — 
the direction of motion of particle j is just an angle 9j, 
calculated from the previous directions of all particles k 
within an interaction range rp = 1 > voAt: 
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where £*■ is a delta-correlated white noise (£ E [— 7T,7r]). 
This introduces a tendency to align with neighboring par- 
ticles, with two simple limits: in the absence of noise, 
interacting particles align perfectly, quickly leading to 
complete orientational order. For maximal noise (rj = 1), 
particles follow random walks. The transition that nec- 
essarily lies in between these two regimes can be charac- 
terized by the following instantaneous order parameter: 
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where N is the total number of particles. 

Varying either the noise strength r\ or the particle den- 
sity p = N/L 2 in periodic domains of linear size L, Vic- 
sek et al. found that (tp) varies continuously across the 
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FIG. 1: Onset of collective motion in cohesion-less models 
Q (original VM, circles) and Q (vectorial noise, squares). 
Variation of order parameter ip (a) and Binder cumulant G 
(b) with the noise strength rj. (vo = 0.5, L — 32, p = 2, and 
equivalent statistics for both models.) 



transition, suggesting the existence of a critical point [8|. 
Studying finite-size effects, they estimated a set of scal- 
ing exponents. Interested in assessing the universality 
of these results and possibly improving these estimates, 
we first introduced simple modifications of the original 
VM such as changing v$ or adding a repulsive force be- 
tween particles to give them a finite extent. Using the 
finite-size scaling Ansatz appropriate for XY-model like 
systems, domain sizes, and particle numbers similar to 
those used in yj, but with much better, well-controlled 
statistics, we were only able to estimate a roughly co- 
herent set of critical exponents after allowing for rather 
strong corrections to scaling (liij . 

For modeling reasons, we also changed the way noise 
is incorporated in the system. In particles make an 
error when trying to take the new direction they have 
perfectly calculated ("angular noise"). One could argue 
that, rather, errors are made when estimating the in- 
teractions, for example because of a noisy environment. 
This leads to change Eq.Q into, e.g.: 
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where n'- is the current number of neighbors of particle 
j. In this case of "vectorial noise", the onset of collective 
motion is discontinuous: for large-enough system sizes, 
(<p) jumps abruptly to zero as rj is decreased, whereas 
it varies smoothly in the original VM (Fig. ^i). This 
is perhaps best seen from the behavior of the so-called 
Binder cumulant G = 1 - (^ 4 }/3((^ 2 ) 2 (Fig.^Ja). In the 
case of vectorial noise, G falls to negative values near r] c , 
the sign of a discontinuous transition, together with the 
phase coexistence expected then. 

Going from angular to vectorial noise is indeed a 
less innocent modification than those mentioned earlier: 
in model J3J), locally-ordered regions are subjected to 
weaker noise than disordered ones. However, it was un- 
clear to us what precisely would be the mechanism to 
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FIG. 2: Discontinuous character of the onset of collective mo- 
tion in the original VM at p = | . (a) : G vs rj at various sys- 
tem sizes, (b): transverse density (bottom curve) and order- 
parameter profile (top curve) in the ordered phase (L = 1024, 
n = 0.18) (c): probability distribution function (PDF) of tp* 
near the transition point, t £ [r;500r], here the correlation 
time is r ~ 10 [3, L = 512. (d) snapshot of coarse-grained 
density field in disordered phase at threshold, p — 2, L = 256. 
The arrows indicate the direction of motion of dense, ordered 
regions. 



change the order of the transition upon introducing this 
nonlinear term. Considering in addition the strong cor- 
rections to scaling found with angular noise, we strived 
to reach larger system sizes in some of these cases, albeit 
at the cost of statistical accuracy 13]. The conclusion of 
these numerical efforts is that the transition is discontin- 
uous in all cases, with finite-size effects being somewhat 
weaker at low densities. As an example, the behavior of 
G with increasing system size shown in Fig. [2Jt for the 
original VM at p = | reveals the characteristic fall to 
negative values. The distribution function of Lp l is bi- 
modal around threshold, without any intermediate uni- 
modal regime (Fig. [5];). Thus, the continuous transition 
reported by Vicsek et al. is only apparent. 

In the ordered phase, the particles are organized in 
density waves moving steadily in a disordered "vapour 
pressure" background of well-defined asymptotic density 
(Fig. • These solitary waves become metastable to a 
long-wavelength longitudinal instability below the den- 
sity threshold rj c (defined to be located at the minimum 
of G), leading to an hysteresis loop. At threshold and be- 
low, the disordered phase consists of nucleated ordered 
patches competing in space and time (Fig.|2Ji). 

At threshold, in the disordered phase, a universal non- 
trivial algebraic scaling law is nevertheless found: the 
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superdiffusive behavior of particles already reported by 
us in 14] is valid in all cases. Trajectories then con- 
sist of "flights" , occuring when a particle is caught in a 
moving ordered patch, separated by normal diffusion in 
the disordered regions. The mean square displacement 
of particles (Sr 2 ) varies like t a with a = 1.65(5). 

We now turn to the onset of collective motion in the 
presence of cohesion. As shown in ^2], the cohesion of 
a population of particles can be maintained without re- 
solving to long-range or global interactions. In the spirit 
of the VM, and following 0, a two-body short-range 
interaction force competing with the alignment tendency 
is introduced, leading to the following model: 
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where a and control the strength of alignment and 
cohesion, 9j k is the direction of the vector linking particle 
j to particle k. The interaction force between these two 
particles, of amplitude /j fe , is actually repulsive up to an 
intermediate equilibrium distance r c , with a short-range 
hard-core at r c and attractive up to the interaction range 
tq. In the following, as in |l2j . we used: 
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if r jk < r c , 
if r c < r jk < r a 
if r a < rjk < r 



(5) 



where rjk is the distance between j and k, with r c 

0. 2. r c = 0.5, and r a = 0.8. Note that vectorial noise 
was chosen in flfl. in the hope of reaching asymptotic 
properties more easily. 

The above model has three main parameters, a, 0, 
and 77, only two of which are independent. The phase 
diagram in the {a, 0) plane (with 77 = 1 fixed arbitrarily) 
was presented in |12j . where, moreover, only neighbors in 
the Voronoi sense are considered in the sums of (@J). For 
large-enough 0, cohesion is maintained, even in the zero- 
density limit. This "gas/liquid" transition is followed, at 
larger values, by the onset of positional (quasi-) order, 

1. e. a "liquid/solid" transition. For large a, these liquid 
or solid cohesive groups move, whereas they remain static 
(up to finite-size fluctuations) for small a. 

In the "liquid case" (intermediate values), the on- 
set of motion is accompanied by a loss of cohesion: while 
small groups set in motion smoothly without breaking up 
(Fig.Ek, dashed lines), larger groups gradually subdivide 
into several parts of roughly equivalent size linked by fila- 
mentary structures, in contrast with their more compact 
shapes before and after onset (Fig. The filaments 
themselves are quite static (Fig. 0Ji) but are displaced by 
the subgroups which move coherently so that they even- 
tually break up, as indicated by the dip in the normalized 
largest connected cluster size n/N in Fig. OJi. Increasing 
a, large groups follow the same precursor of the tran- 
sition as smaller groups, but when their fragmentation 
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FIG. 3: Onset of motion of cohesive groups in model Q 
with 77 = 1, Vq = 0.05. (a) (ip) and n/N (normalized size 
of largest connected cluster) vs a (p = (3 = 20 (liquid 
phase), dashed lines: N = 4096; solid lines: N = 16384). In- 
set: solid group (f3 — 84) of N = 4096 particles; dashed line: 
(ip); solid line: relative diffusion of initialy neighboring parti- 
E fc ~ 3 (l - r%(t)/r%(t + T))) jlt where T « 20iV 
1 in the liquid phase, while A ~ in the solid phase, 
see |l3])- (b) variation with a of the maximal absolute ro- 
tation angle \uo\ averaged over 100 samples of 1000 vortices 
(N = 2048, p = fJ = 30 (liquid phase)). Dashed line:(v3> 
during the same runs. Inset: distribution of rotation times at 
the transition with decay exponent ~ 1.3. 



occurs the order parameter falls back, leaving an inter- 
mediate peak (around a = 1.7 in Fig. |3]i). Increasing a 
further, (ip) rises again and finally jumps to (ip) — 1 when 
full cohesion is recovered (for a = 1.88 in Fig. |3Ji). This 
discontinuous jump is the true location of the transition: 
For an infinite group, the onset of motion must occur 
abruptly near this value, as the precursory features de- 
scribed above disappear because the population divides 
into infinitely-many subgroups whose influences average 
themselves out. Meanwhile, cohesion is only lost at the 
transition point in this asymptotic picture. 

The breakup of large cohesive groups around threshold 
is probably closely related to what happens in the case 
without cohesion: the subgroups connected by filaments 
may correspond to the ordered patches seen in the dis- 
ordered phase near threshold in Fig. The breakup 
itself can be seen as resulting from the maximal effect 
of acoustic modes on the shape of the group 01 ■ Also 
affecting the shape dynamics are rotational modes: the 
subgroups seen in Fig. [2^, not only move but they also 
rotate slowly |l6| . Rotation is not steady, but intermit- 
tent. We recorded the rotation times and their corre- 
sponding angles. Extremal statistics analysis reveal that 
the tendency to rotate is maximal at the onset of motion 
(Fig. [3t>) • Moreover, at threshold, the distribution of ro- 
tation times is algebraic with a decay exponent such that 
it has no finite mean (inset of Fig. Et>) ■ 

The onset of motion of the "solid" groups (large val- 
ues) is accompanied by a loss of positional order: these 
crystals melt near the transition (inset of Fig.|3^). Given 
the above results in the liquid case, one can expect very 
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assays" consisting of grafted molecular motors such as 
kinesin (resp. myosin) moving filaments made of tubulin 
(resp. actin) might provide the simplest setting in which 
to investigate superdiffusion at onset, given the available 
observation techniques |a, |6j . 
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FIG. 4: Typical shape of a liquid cohesive group of 16384 
particles (model Q, p = /3 = 20 arrows indicate direction 
of motion), (a): at onset before loss of cohesion, a = 1.78. 
(b): static phase, round shape, a = 0.5. (c): in moving phase, 
typical triangular form (see [12]), a = 2.5. (d): close-up of a 
filament: no local order is apparent. 



large solid groups to melt and then subdivide and lose 
cohesion in the transition region, making the onset of 
motion asymptotically discontinuous. 

To summarize, the onset of collective motion in the VM 
as well as in related models with and without cohesion is 
always discontinuous, and the critical behavior reported 
in Q is only apparent and due to (strong) finite-size ef- 
fects. Without cohesion, the ordered phase consists in 
density waves propagating steadily in a disordered back- 
ground. With a short-range repulsion/ attraction inter- 
action, the cohesion ensured both in the disordered and 
ordered phases is broken at the onset of motion under the 
competing influence of sound, density, and shape modes. 
The resulting mesoscopic subgroups are linked by fila- 
ments which, however, we believe to be probably non- 
universal, model-dependent structures. 

At the theoretical level, ongoing work is directed to- 
wards the understanding of the complex interplay be- 
tween shape (surface tension) and acoustic modes, and 
of the stability properties of the density waves. At the 
experimental level, it remains difficult to study quanti- 
tatively bird flocks and fish schools, and moreover we 
have no specific prediction as to the onset of motion of 
these cohesive groups 0. Without cohesion, however, 
the universal superdiffusive behavior observed in the dis- 
ordered phase near threshold could be observed experi- 
mentally. As already suggested in bacteria such as 
E. Coli might be good self-propelled particles. Human 
melanocytes also look promising in this respect as shown 
remarkably by the group of Gruler Q . Finally, "motility 
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In the other hand, interesting predictions were made by 
Tu and Toner |9j] for the ordered phase. 



